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Optical Feshbach resonances (OFRs) have generated significant experimental interest in recent 
years. These resonances are promising for many-body physics experiments, yet the practical appli¬ 
cation of OFRs has been limited. The theory of OFRs has been based on an approximate model 
that fails in important detuning regimes, and the incomplete theoretical understanding of this effect 
has hindered OFR experiments. We present the most complete theoretical treatment of OFRs to 
date, demonstrating important characteristics that must be considered in OFR experiments and 
comparing OFRs to the well-studied case of magnetic Feshbach resonances. We also present a 
comprehensive treatment of the approximate OFR model, including a study of the range of va¬ 
lidity for this model. Finally, we derive experimentally useful expressions that can be applied to 
real experimental data to extract important information about the resonance structure of colliding 
atoms. 


I. INTRODUCTION 
A. Background 


Magnetic Feshbach resonances (MFRs) have become 
a staple of quantum gas experiments with alkali-metal 
atoms, allowing for unprecedented control of interatomic 
interactions [lj. The MFR technique is so powerful that 
it has extended the reach of dilute quantum gas exper¬ 
iments to a variety of areas of physics. Examples of 
high impact experiments that utilize MFRs are the study 
of strongly correlated systems 0] such as unitary Bose 


;ases, the discovery of exotic few- 


01 and Fermi 
body bound states 
molecules 12, 13], and the engineering of novel quan¬ 
tum matter |14, JL51. Feshbach resonances based on laser 


llj, the ability to make ultracold 


fields^ known as “optical Feshbach resonances” (OFRs) 
[III 03 — have also been observed [l8|,[l9[, but so far their 
utility has been limited. Since laser fields can be focused 
tighter and switched faster than magnetic fields, it is ex¬ 
pected that OFRs could yield an MFR-like effect but 
with orders of magnitude better spatial and temporal 
control El. 


Furthermore, OFRs are better suited for alkaline- 
earth-metal atoms, which have magnetically insensitive 
electronic ground states. The study of alkaline earth 
atoms is now a rich field, attrac ting attention for metrol¬ 
ogy f2ll ] , quantum information [22L [23[, and many-body 
physics 04|. Quantum degenerate gases of these atoms 
have also been realized |25-27]. Many-body physics has 
been demonstrated in strontium lattice clocks |28l |. and 
it has been shown that controlling many-body interac¬ 
tions in gases of alkaline earth metals could lead to better 
clock accuracy 09] . Without MFRs to facilitate the same 
many-body control enjoyed by alkali-metal experiments, 
OFRs have been suggested as an alternative |30| . 


OFRs have been the focus of several experiments. 
These resonances have been observed in alkali gases 
[T1 '31] and in alkaline-earth-like atoms [32]. P-wave 
OFRs have been reported 03] , and OFRs have been suc¬ 
cessfully applied to induce thermalization in Sr gases 0H 
and manipulate the condensate dynamics of a Sr Bose- 
Einstein condensate [35, 36]]. The theory used to describe 
these experiments was based on a quantum defect treat¬ 
ment by Bohn and Julienne, who used an isolated res¬ 
onance approximation to derive the optically modified 
scattering length 07]. 


Although the isolated resonance theory has been suc¬ 
cessful in describing some observations of OFRs, it fails 
to explain OFR behavior in the large detuning regimes 
that are critical to a proposal for practically applying 
these resonances 0fj . Attempts to exp erimentally re¬ 
alize this proposal did not succeed 03; therefore, the 
limited theoretical understanding of OFRs has hindered 
experimental progress. To broaden the theoretical under¬ 
standing of OFRs, we perform the most complete theo¬ 
retical analysis of this effect to date. To this end, we treat 
OFRs with a numerical coupled channel method, which 
has been highly successful for treating MFRs [lll38|. Like 
the MFR theory, our coupled channel approach is capa¬ 
ble of treating multiple interacting resonances without 
being restricted to the more limited isolated resonance 
approximation. Consequently, this more general treat¬ 
ment allows us to study the range of validity of the iso¬ 
lated resonance approximation, and it also enables us to 
point out similarities and significant differences between 
OFRs and MFRs. Finally, we use the isolated resonance 
theory to derive experimentally useful formulas that can 
be used to understand real experimental OFR data. 
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B. Basic Collision Theory 

In the context of cold-atom physics, a Feshbach reso¬ 
nance is a collisional resonance of two particles that is 
tunable by an external field. This is possible if a molecu¬ 
lar bound state from an excited scattering channel (called 
the “closed channel”) couples to the free atom contin¬ 
uum of the ground state scattering channel (called the 
“entrance channel” or the “background channel”). Fur¬ 
thermore, the bound state energy is tunable by an ex¬ 
ternal magnetic or electromagnetic field. The presence 
of this bound molecular state modifies the s-wave scat¬ 
tering length of the atoms, thereby changing the atomic 
interactions as the external field is tuned. 

We emphasize that both MFRs and OFRs can be 
treated by the same scattering formalism, which accounts 
for the differences in their coupling and control mecha¬ 
nisms, as presented in the review by Chin et al. [0] . A 
typical MFR is coupled to the entrance channel by in¬ 
ternal short range spin-dependent couplings within the 
ground state manifold of Zeeman levels. MFRs are tuned 
by varying an external magnetic field B to move a molec¬ 
ular bound state across a collision threshold, creating a 
pole in the scattering S'-matrix as a function of B. An 
OFR involves coupling a bound molecular state to two 
colliding atoms in their ground states using a photon 
from a laser, hereafter called the “photoassociation laser” 
or “PA laser.” In this case, the coupling depends on both 
the PA laser detuning from a photoassociation resonance 
(the “molecular detuning”) and the PA laser intensity. 
In contrast to MFRs, which are often based on molec¬ 
ular states that have very long lifetimes, spontaneous 
decay of the excited molecular state in an OFR intro¬ 
duces an appreciable linewidth to the molecular transi¬ 
tion. Any population transferred to the excited state un¬ 
dergoes spontaneous decay, which translates to inelastic 
loss collisions that must be minimized in order to utilize 
an OFR. However, resonance decay does not necessarily 
prevent the application of OFRs since MFRs with 2-body 
decay channels [391441]] have proven experimentally use¬ 
ful 

In the OFR studies presented here, we consider bosonic 
88 Sr with the two ground state atoms providing the 
+ 1 So ground state entrance channel and the excited 
state 1 S'o + 3 Pi providing the closed channels, schemat¬ 
ically represented in Fig. [T] Since bosonic isotopes of 
alkaline earths have no nuclear spin, the 88 Sr resonance 
structure is considerably simpler than for atoms with hy- 
perfine interactions, making it a good atom for an OFR 
experiment. The — > 3 P\ atomic transition is an inter¬ 

combination line with a natural linewidth of y a = 27rx 
7.4 kHz. The narrowness of this transition means that 
all OFRs in 88 Sr are well resolved from the atomic line, 
decreasing the severity of off-resonant atomic light scat¬ 
tering. 

To analyze the scattering of two colliding 88 Sr atoms 
in a light field, we calculate the scattering S'-matrix to 
determine the elastic and inelastic scattering cross sec¬ 


tions. Since the 88 Sr ground state is completely spinless, 
and since current OFR experiments are typically per¬ 
formed at temperatures of a few g K or below, the scat¬ 
tering is described by an s-wave collision with a single 
nondegenerate entrance channel. Therefore, we will de¬ 
velop our theory for this experimentally simple OFR sys¬ 
tem, for which we only need a single s-wave S'-matrix ele¬ 
ment S(k) = e 2lr, ( k \ represented in general by a complex 
energy-dependent phase shift rj(k). Here k is defined via 
the collision velocity Hk/fi , n = m/2 is the reduced mass, 
and m is the mass of an 88 Sr atom. This phase shift in 
turn defines an energy-dependent scattering length a(k) 

as mmm, 


(U\ (u\ -utu\ tan? ? (/c) 1 1 - S(fc) 

<.(*) = «(*)-.**) = — — = a) 

This expression is useful for small but nonvanishing colli¬ 
sion energies, and the standard complex scattering length 
is the k —> 0 limit of this expression. The elastic and in¬ 
elastic loss cross sections are 

a el = ^\l-S(k)\ 2 = 8n\a(k)\ 2 f 2 (k), (2) 

a in = ^(l-\S(k)\ 2 )=^b(k)f(k). (3) 

Here g is a collisional symmetry factor, which is equal to 
2 for identical bosons (as assumed here). The function 


1 + k 2 \a(k)\ 2 + 2kb(k) 


approaches unity when k\a\ -C 1 for all detunings. For a 
trapped gas of atoms, this limit occurs when ksT/h'y -C 
1, where ks is Boltzmann’s constant and T is the sam¬ 
ple temperature. The elastic and inelastic collision rate 
coefficients are related to these cross sections as 

hk h 

K e i(k) = — a ei (k) -A 8tt— k\a(k)\ 2 as k —> 0 (5) 

k h 

Hk h 

K in (k) = — ain{k) —> 8n—b(k) as k —> 0 . (6) 

h h 

These general expressions are valid in the s-wave limit for 
OFRs and for decaying or non-decaying MFRs. A sum 
over higher partial waves is needed when these begin to 
contribute at higher k, and a thermal average of K e i(k ) 
and K in (k) is needed when comparing to experiment. 


II. COUPLED CHANNELS FORMULATION OF 
OPTICAL FESHBACH RESONANCES 

A. Background 

The standard treatment for atomic collisions involving 
two or more internal states of the atoms is the coupled 
channels (CC) method. Numerical models based on CC 
methods have been very successful in treating collisions 





3 



FIG. 1: (Color online) a) The 1 So + x So entrance channel of 
88 Sr couples to a bound state of the 3 Pi + 1 So closed channel 
via the PA laser field. The atomic transition is a 7.5 kHz 
intercombination line. Here E is the collision energy, u>o is 
the atomic resonance frequency, and uj is the laser frequency, 
b) In the dressed state picture, two free atoms in the entrance 
channel are brought into resonance with a molecular bound 
state. Here S is the “molecular detuning” (defined in Section 
DU]). and n p is the photon number. The Condon radius R c is 
defined as the value of R where the two potentials cross. 

and MFRs of ground state alkali-metal atoms PHHI]. The 
CC method involves setting up a basis set representing 
the “channels” of the electronic, spin, and rotation de¬ 
grees of freedom of the colliding atoms and then solving 
the matrix Schrodinger equation for the amplitude of the 
radial motion in the interatomic separation coordinate R 
for each of these channels. 

In MFR theory, the channels represent the states of 
the atoms in a magnetic field for R —> oo, the Born- 
Oppenheimer potentials characterize the P-dependent 
interactions, and spin coupling matrix elements are ap¬ 
proximated by their atomic values. An external B field 
shifts the energies of the atomic and molecular energy 
levels. In the case of OFRs, the channels represent 
the field dressed atoms, where the ground and excited 
states are coupled by the light field in the dipole approx¬ 
imation, and the P-dependent molecular interactions 
are represented by the ground and excited state Born- 
Oppenheimer potentials together with any non-adiabatic 
coupling between them. The OFR case has the added 
complication of including the spontaneous emission of 
light by excited state atoms or molecules. 

To date, all cold atom OFRs have been treated by a 
resonant scattering formulation [l6|, |37j, EH . The next 
section will discuss the approximation of treating each 
OFR as an isolated resonance. Here we concentrate on 
giving a full CC treatment f50l - [53j that includes the ef¬ 
fect of multiple overlapping resonances without restrict¬ 
ing the theory to treating isolated single resonances. This 
enables us to establish the conditions under which the 
isolated resonance approximation is valid. 

We follow the Held-dressed collision approach of Juli¬ 
enne [mI . EH , which was applied to explain experiments 
on the collisional redistribution of light EE E3|. To 


do this it is necessary to properly represent the three- 
dimensional (3D) nature of the_collision and the role of 
atomic degeneracy. References EE I 58)] treat the exchange 
of multiple photons during a collision, by which one par¬ 
tial wave is coupled to higher partial waves through the 
intrinsically anisotropic nature of the interaction with 
light. Reference EE adapts the CC dressed atom for¬ 
malism to cold atom collisions in strong optical fields to 
explain the phenomena of optical shielding. 

Three effects need to be incorporated within a CC the¬ 
ory to describe OFRs in the collision of 1 So Sr atoms in 
a light field tuned near the 1 Sq — > 3 P± line: (1) the field 
dressing of the collision, (2) the inherently 3D nature of 
the collision, with a space axis defined for the separated 
atoms by the PA laser polarization but with a rotating 
interatomic axis needed for the excited molecular bound 
states, and (3) the spontaneous emission while in the ex¬ 
cited state. If we make the approximation that the light 
field is weak, the total angular momentum J is a good 
quantum number (the optical coupling matrix element re¬ 
mains small compared to the spacing of rotational levels 
in the excited state). In this case Refs. EE [57] showed 
that six CCs are needed to describe optically coupled 
1 So+ 1 So —» 1 So+ 1 Pi collisions. The same is true when 
we replace 3 Pi with 3 Pi. One set of channels represents 
the ground state collision with partial wave £ = J and n p 
photons at an angular frequency w. Another set repre¬ 
sents the excited 0„ and l u molecular states with n p — 1 
photons at frequency oj and respective projection S2 = 
0 and 1 of electronic angular momentum j = 1 on the 
interatomic axis. These excited state channels have total 
angular momentum J e = J — 1 (two channels), J (one 
channel), and J +1 (two channels). In the special case of 
s-wave collisions (J = 0) of cold atoms, only three chan¬ 
nels are needed, representing the ground state and the 0 U 
and l u states with J = 1. Finally, spontaneous emission 
from the excited molecular state can be included with a 
complex potential EH , with a caveat that the imaginary 
decay part of the potential has to be turned off when the 
atoms are far apart in the free atom limit. 

We assume that the free atoms are weakly dressed— 
that is, the PA laser with photon energy fiw is detuned 
from the atomic excitation energy hu>o by a large amount 
compared to the optical coupling strength 

V opt = (2irI/c) 1/2 d, (7) 

where I is PA laser intensity, c the speed of light, and d is 
a molecular transition dipole matrix element [57J. How¬ 
ever, the short range molecular states can be strongly 
dressed, so that the peak of an on-resonant PA line at 
hu) n , where n is the molecular vibrational level, can be 
power broadened. The rotational quantum number J e 
will remain a good approximate quantum number as long 
as V),pt remains small compared to the rotational con¬ 
stant B n of level n. (The separations of the J = 0 and 2 
levels from the J = 1 level are 2 B n and 4 B n respectively.) 
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TABLE I: Minimal CC basis sets in the Hund’s case b = 
(c) and (e) representations, where a = 0, ±1 represents the 
polarization of the light with n p photons of frequency uj that 
couples the ground and excited states. 


Channel Case (c): \Cl g JM)\n p oJo 

-} Case (e): \jtJ M)\n p u>a) 

|1>6 

0 ff 00) c \n p uj(j) 

0 9 000) e \n p uja) 

|2 )b 

0„1 a) c \n p — l,cja) 

|l„01cr) e |n p — l,cua ) 

]3>b 

|l«lcr)c| n p — lc va) 

l u 21o-) e |n p - Icoa) 


B. Formulation for 88 Sr 


We include in our treatment here the minimal num¬ 
ber of three channels needed to get a basic description 
of near-threshold s-wave OFRs. This minimal treatment 
could be written in either of two basis sets representing 
the electronic, spin, rotational, and photon degrees of 
freedom. One basis set for the molecular degrees of free¬ 
dom is the Hund’s case (c) basis represented as |f l s JM) c , 
where the projection of electronic plus spin angular mo¬ 
mentum j on the rotating body-fixed axis is Cl, s rep¬ 
resents the gerade or ungerade inversion symmetry of 
electronic coordinates, and M is the projection of total 
angular momentum J on a space-fixed axis. The other 
molecular basis is the Hund’s case (e) asymptotic basis of 
Refs. @,[13 represented as \j s lJM) e , where j s = 0 or 1 
represents the separated atoms in the respective 1 S'o+ 1 5'o 
and 1 S , 0 + 3 Pi channels with partial wave i, coupled to 
total angular momentum J and projection M. The sub¬ 
script s on j indicates that the electronic wavefunction is 
symmetrized with respect to the exchange of electronic 
coordinates. Table [Q shows the three basis functions for 
a dressed CC calculation in either representation. The 
transformation between the molecular and asymptotic 
representations is (see, for example, Eq. (36) of Ref. [H|) : 

|0 „JM) C = (——J \l u ,J-l,JM) e 

/ j + l \ 1/2 

~U/Tlj ( 8 ) 

/ J + 1\ 1/2 

|1 u JM)c = ( 2 J _|_ 1 ) |lu> J ~ 1 ) JM)e 

( J \ 1/2 

+ ^ j _|_ i J Il«,i+l,JM) e . (9) 


Using the CC expansion of the full wavefunction at 
total energy E, 


3 

*(R, E) = J2 I i)bF iib (R, E)/R (10) 

i= 1 

where E)^ represents the amplitude of the wavefunction 
projected on the basis function |*)f, for Hund’s case b = 
(c) or (e). The CC matrix Schrodinger equation for the 


s-wave collision of the two atoms in a (moderately) weak 
light field is 


d 2 T 

9R2 + £2 


(£M-V(R))tf = 0 


( 11 ) 


where I is the identity matrix and the potential matrix V 
describes the diagonal and off-diagonal matrix elements 
of the collisional and optical interactions. Either the b = 
(c) or (e) representations (Table Q]) of the excited state 
could be used to set up the expansion and V matrix in 
Eqs. (gUP and (fill) . Our numerical calculations use the 
Hund’s case (e) basis, for which the matrix elements are 


given in Table I of Ref. [57J, and quoted in the supple¬ 


mental online material for Ref. 


/ 


V = 


v a 

v opt 


V o 


V r 


opt 


l(V 0u + 2V lu ) Y(V lu -V 0u ) 
*§(V lu -V 0u ) h^V 0u + V lu )+6V ct 


_ ( 12 ) 

where the 614en term represents the d -wave centrifugal 
potential with V cen = h 2 /(2^R 2 ). Here V g (R), Vq u (R), 
and Viu(R) represent the ground state and and 1„ 
excited state BO potentials, each of which we model as 
a Lennard-Jones potential plus an additional long range 
term: 


V g (R) = 



goo 


(13) 


+ VU14) 


+ Ko{15) 


The V soo terms give the asymptotic values of the poten¬ 
tials as R —> oo, as explained below. We use the excited 
state potential parameters from Zelevinsky et al. [13, Hi] ■ 
The ground state Ce i9 and Cs j9 parameters come from 
Ref. [62|, and i?o,g w &s optimized to reproduce the mea¬ 
sured bound state binding energies of Ref. S63] to within 
0.4% 64iJ. V g (R) has an s-wave scattering length of —1.4 
ao, consistent with that reported in Ref. [631 ]. 

The optical coupling matrix element in Eq. m is 
given by Eq. o in the dipole approximation, where we 
neglect retardation (that is, the phase difference between 
the optical fields separated by distance R <C A, where 
A = 27 tc/lo is the wavelength of the excitation light). 
Thus, since we use the symmetrized g and u electronic 
states, d = \f2d a , where the atomic transition dipole 
d a = 0.08682 atomic units (1 a.u. = eao, where e is the 
electron charge and ao is the Bohr radius), corresponding 
to an atomic 3 Pi lifetime of 21.46 fis or linewidth of j a 
= 27rx 7.416 kHz. Thus, introducing units into Eq. 0, 

V 0 pt/h = 24.83MHz x d a \Jl/{ lW/cm 2 ), (16) 

where d a is in atomic units. The optical coupling in 
V conforms to the case (e) selection rule that A£ = 
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0, A me = 0 (it is only the electronic j quantum num¬ 
ber that changes). This coupling is also independent of 
light polarization a for this transition. Using Eq. m, 
there will be an asymptotic light shift 

(IT, 

which is negative for the ground state and positive for 
the excited state. Thus, taking V^oo = Vx, and V uoo = 
h(uJo — uj) + Voo in Eqs. (fl3|) - (|T5]l ensures that when V 
is diagonalized the lowest energy eigenvalue at large R 
for the field-dressed ground state is zero. With this def¬ 
inition of the zero of energy, the total energy E in the 
CC Schrodinger equation m represents the relative col¬ 
lision kinetic energy h 2 k 2 /2p, of the dressed ground state 
atoms, and E = h 2 k 2 /2g —> 0 at the collision threshold. 

The matrix V in Eq. (fl2l) could be transformed to the 
molecular case (c) representation by transforming the 2 x 
2 excited state block using the (c) to (e) transformation 
matrix used in Eqs. m and ©. This would give the 
diagonal J = 1 0„ and l u potentials given in Eqs. (1) 
and (2) of Zelevinsky et al. 0 and generate the body- 
frame Coriolis coupling term between these two states. 
The optical coupling in the case (c) molecular basis is 
different from that in the asymptotic case (e) basis. Using 
the transformations in Eqs. 0 and © shows that the 
optical couplings matrix elements between the ground 
J = 0 0 g state and the respective excited J = 1 0 U and 
l u states are determined from Eq. 0 with the molecular 
dipole matrix elements 

d 0u = y/l/3V2d a , (18) 

d lu = ^/2j3V2d a . (19) 

The \[2 is the same homonuclear g u enhancement 
factor that affects the case (e) matrix element. The other 
factor corresponds to the usual Honl-London factor for 
R-branch (J —> J + 1) molecular transitions. 

Treating an OFR requires that we include the decay 
from the excited state, which has a molecular decay rate 
7 . Our calculations assume 7 = 7 m , where we define 
7 m = 27 q . This rate 7 m is the rate of spontaneous emis¬ 
sion from the excited molecular state in the long-range 
non-retarded dipole approximation. A nonzero value of 
7 _ 7 m would be due to other processes that induce de¬ 
cay of the excited state or change the emission rate from 
its long-range non-retarded dipolar value. While Bohn 
and Julienne [uP introduced artificial channels to sim¬ 
ulate excited state decay, a simpler method is to intro¬ 
duce an imaginary term in the excited state potentials in 
Eqs. (fl4ll and (fl5l) . replacing Vj u , j = 0,1, with: 

V JU ~ i y (l + e^ R ~ R ^ _1 , (20) 

where 7 is an i?-independent constant. The function 
in parenthesis ensures that molecular decay turns off at 


large distances when R exceeds the arbitrary cutoff ra¬ 
dius Rcut by an amount large compared to the length 
1 /(3. Furthermore, this function ensures that the full 
molecular decay rate turns on at small distances where 
Rcut. — R is appreciably less than 1//3. The constant /3 
parametrizes the distance over which molecular decay be¬ 
comes appreciable. 

When the coupled equations are solved with this com¬ 
plex potential in Eqs. (m and (USD, the 5-matrix is 
non-unitary, and 1 — |5(fc)| 2 in Eq. 0 represents loss of 
ground state atoms due to molecular excitation followed 
by excited state decay. We assume that every sponta¬ 
neous emission event represented by the imaginary term 
in Vj U results in hot atom or molecular products that are 
lost from the trap. Our numerical studies show that this 
assumption is good for all the excited levels except the 
state nearest in energy to the atomic resonance (Section 
Iff Dll . We calculate that 60% of the emission from this 
state does not result in loss from a 10 pK trap [60l l65j. 
The cutoff ensures that there is no spurious excited state 
decay associated with the asymptotically dressed atoms. 
We find that in the core of a PA line, out to molecular 
detunings of several hundred line widths from molecular 
resonance, the loss associated with the imaginary part of 
the scattering length is not sensitive to the value chosen 
for Rcut, as long as it is significantly outside the outer 
turning point of classical motion for the excited state vi¬ 
brational level. Furthermore, the real part of the scatter¬ 
ing length is completely insensitive to the choice of i?, cu t ■ 
We typically choose J? C ut = 500 ao and /3 = 0.05 a f y . We 
find that our numerical calculations were insensitive to 
the choice of /3. 


C. Approximations and limitations 

This three-channel model makes several approxima¬ 
tions, but is capable of representing the essential qualita¬ 
tive and semi-quantitative effects associated with OFRs 
in the weak to moderate held regime. Our model only 
includes the minimal number of partial waves needed to 
represent the change in scattering length and molecular 
losses due to the OFR. ft leaves out the coupling of the 
excited J = 1 levels to the ground state d -waves as well 
as coupling to the doubly excited states associated with 
the 3 Pi + 3 Pi separated atom limit. This means that the 
light shifts calculated from the three-channel model will 
not be accurate, since ground state d-waves are known 
to contribute to PA light shifts 66. 67], and doubly ex¬ 
cited states may contribute also. Furthermore, the effect 
of field-dressing in modifying the ground state thresh¬ 
old elastic scattering of partial waves with i > 0 is not 
included. This modification is due to field dressing that 
brings in 1 /R 3 terms in the long range potential that will 
affect the d-wave collisions of like bosons or the p -wave 
collisions of like fermions or unlike species. Note that 
in our three-channel treatment, the field dressed s-wave 
interactions have the correct property that they have no 
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long-range 1 /i ? 3 component, since such variation van¬ 
ishes in the V 2.2 matrix element of Eq. (fl2l) (if we had 
attempted only a two-channel field dressed treatment, 
the presence of the single 0 U excited state potential would 
have introduced a spurious 1 /i ? 3 term in the ground state 
dressed s-wave potential). 

It would be straightforward to introduce higher par¬ 
tial waves and strong field dressing into the calculation, 
using the formalism of Napolitano et al. [58j. This for¬ 
malism uses the “uncoupled” asymptotic basis \ jrrijime), 
which is better for treating strong field dressing than 
the “coupled” | jlJM) basis we use here. The subtle 
effects of retardation, switching off the dipole approxi¬ 
mation, and including the weak coupling to the gerade 
states as R increases [fill should be taken into account 
in a more complete theory. We do not perform a time- 
domain analysis, so we cannot reproduce the transient 
OFR dynamics jfil] observed in Ref. [fifij]. Furthermore, 
a full treatment of excited state spontaneous emission 
during a collision is beyond the scope of CC methods, 
and would require treatment by stochastic Schrodinger 
equation methods (density matrix methods are compu¬ 
tationally intractable) |7(j, [nj. Last, as we show in the 
next section, our analysis with a cutoff of the long range 
decay is sufficient for treating OFRs for molecular detun¬ 
ings that are of order 100 line widths (or less) from the 
center of a PA line. 



FIG. 2: (Color online) Real (solid line) and imaginary (dashed 
line) parts of the complex scattering length a = a — ib. Here 
the detuning from atomic resonance u—vo (the “atomic detun¬ 
ing”) is measured with respect to the 88 Sr intercombination 
line transition at vq. Also, E/ks = 4 gK and I = 10 W/cm 2 . 
The dotted line shows the background dbg for E/kB = 4 gK 
and 1 = 0. Inset: A close up of the off-resonant behavior of 
the n = —2 OFR (solid line). Also plotted is the scattering 
length predicted by treating the n = —2 OFR as an isolated 
resonance (dashed line). 


D. Coupled Channels Results 


The CC calculations to solve Eq. HID were carried out 
using the standard renormalized Numerov method 72} 
using complex variables so as to represent the effect of the 
complex potential in the excited state channels. A step¬ 
doubling algorithm was employed to optimize the number 
of steps needed as R increases between the short and long 
range regions. The single S'-matrix element S(E,I,u) 
for the dressed ground state s-wave was extracted from 
the log derivative of the single open channel solution 
Fi te (R) of the three channel propagated wavefunction of 
Eq. (fTUl) at a suitable large asymptotic value of R. Using 
Eq. JT]) then gives the complex energy-dependent scat¬ 
tering length a(/c,J, w), which then gives the elastic and 
inelastic rate coefficients K e i and K in (Eqns. [5] and [6]). 

Fig. [2] shows the real and imaginary parts of a(k) as 
the PA laser frequency v = w/27r is detuned from atomic 
resonance at vq = ujo/2tt (where v — vq is the “atomic 
detuning”). This particular example was taken for a 
PA laser intensity of 10 W/cm 2 and a relative collision 
kinetic energy of E/ks = 4 /xK. Here the background 
dbg = 0.495 ao at E/ks =4 g,K differs from the back¬ 
ground value -1.4 ao in the limit of E = 0 due to the 
energy dependence of Ob g (fc)- The calculations were car¬ 
ried out for atomic detunings larger in magnitude than 
-20 MHz to avoid strong field-dressing effects at atomic 
detunings near resonance (the optical coupling matrix 
element V/ P t /h = 0.84 MHz for this I). The decay rate 


was taken to be 7 = 7 m = 2^ a = 27tx 0.014833 MHz. 
The figure shows a series of four OFRs in this region. 
These four resonances correspond to the previously ob¬ 
served fH n = -2, -3, and -4 members of the 0„ J = 1 
series at binding energies E n /h = 24 MHz, 222 MHz, and 
1084 MHz and a single n = -1 member of the 1„ J = 1 
series at 353 MHz. Here n < 0 counts bound states down 
from the last level (of a given 0 U or l u symmetry) des¬ 
ignated as n = — 1. The scattering lengths show a series 
of overlapping resonances that cause a large change in 
scattering length near the poles of the resonances but 
return to Ob g between resonances. The fact that the 
stronger n = —2 resonance returns to its background 
value near the n = —3 line (Fig. [5] inset) illustrates an 
important general feature of a vibrational sequence of 
OFRs: interfering resonances cause the scattering length 
to return to its background value in between resonances. 
Even the presence of a neighboring OFR that is compar¬ 
atively weak will diminish the off-resonant magnitude of 
a stronger OFR (Fig. [2] inset). This property imposes a 
constraint on OFR experiments, namely that molecular 
detunings cannot be so large as to be comparable to the 
frequency separation between the resonance of interest 
and the nearest resonance. In contrast, MFRs arising 
from neighboring spin-channel resonances interfere with 
one another in a manner that is qualitatively different 
than for a vibrational series (see Section IIVBI) . 

The imaginary part b of the scattering length that gives 
the loss rate coefficient, Eq. ©, shows a series of sharp 
spikes near the poles of the resonances in Fig. [2j and 
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FIG. 3: (Color online) Imaginary part of a = a — ib from 
Fig- EJshown on a log scale. The diamonds show a 1 /(u — uo) 2 
scaling. The inset shows an expanded view of the 0 U n = —4 
resonance near —1084 MHz. The dashed, solid, and dotted 
lines show the results for R cut = 200 ao, 500 ao, and 1000 ao 
respectively. Near the peak of the resonance, b is independent 
Of Rcut ■ 


shows a value very near zero on the linear scale of the 
figure. Figure [3] provides a better way to illustrate the 
basic features of b by showing it on a log scale. Here the 
“background” on which the poles sit varies as l/(i/— ^o) 2 , 
with this functional form indicated by the diamonds on 
the figure. Furthermore, far detuned from a molecular 
resonance, the magnitude of this background is found to 
scale linearly with R cut as R cu t increases. This is be¬ 
cause away from resonance, most of the loss of flux in 
the collision due to the presence of a complex potential 
comes from the long range region, where decay should 
not be counted as loss, since it merely represents atomic 
light scattering that returns an atom to its ground state. 
Thus, loss is overcounted by use of a complex potential 
if i? C ut is too large. 

In the core of the line spanning molecular detunings of 
1007 m , b and Kin are independent of R cu t- For example 
at a molecular detuning of lOOym (inset to Fig. [3]) the 
values of b calculated with R cnt = 200 ao or 500 ao differ 
by less than %. The difference grows to 10% if R cut = 
1000 ao- The difference with R cut = 200 ao or 500 ao 
only grows to 10% when the molecular detuning is over 
250 line widths. Consequently, if i? cut is selected to have 
a small enough “physical” value where spontaneous de¬ 
cay for R < R cut represents actual loss of atoms, the loss 
calculated for molecular detunings up to a few hundred 
line widths are meaningful and not sensitive to the choice 
of i? C ut- In any case, the scattering length a given by the 
real part of a is completely insensitive at all detunings to 
the choice of R cu t- Comparing our CC theory to exper¬ 
imental data taken at small molecular detunings, we are 
able to reproduce the resonance strengths measured in 
Refs. ;34j and [35]. However, the atom loss rate of Ref. 


'doi is measured at a molecular detuning large enough for 
our theory to be sensitive to R c u t; therefore, our theory 
is not designed to reproduce this rate. 

Within the inherent limitations of our approximations 
that we have outlined above, we expect our CC calcula¬ 
tions to give the correct change in scattering length for all 
detunings and the atom loss rate coefficient for at least 
100 line widths from the peak of a molecular resonance. 
Consequently, since the resonances are spaced by many 
thousands of line widths apart, we turn our attention in 
the next section to understanding the theory for single 
isolated OFRs for molecular detunings in the vicinity of 
a photoassociation resonance. 


III. ISOLATED RESONANCE THEORY OF 
OPTICAL FESHBACH RESONANCES 


The OFR features in Figs. [2] and [3] tend to be well- 
isolated from one another and thus can be described quite 
successfully by theory designed to treat an isolated sin¬ 
gle resonance situated on a background. Isolated reso¬ 
nance theory has been used for cold atom OFRs since 
they were first proposed [HI, [Tt]. [371 IH1 H3 ■ This the¬ 
ory successfully explained alkali-metal atom PA spectra 
with hyperfine structure (in good agreement with experi¬ 
ment [59, 73]), and it also explained the saturation of PA 
lines at high intensity { 6 (| [74j, |H1- In fact, both OFRs 
and MFRs can be treated by the same resonance scatter¬ 
ing formalism when the possibility of decay of the closed 
channel resonance state is taken into account |lj. 

The isolated resonance approximation assumes that 
each molecular bound state is far from the other molecu¬ 
lar states in the closed channel and can be described by a 
strength parameter that is local to the resonance—that 
is, independent of energy and molecular detuning in the 
vicinity of the resonance. Bohn and Julienne give a gen¬ 
eral resonance scattering treatment for an OFR based 
on quantum defect theory (37J. They derive a general 
expression for the S'-matrix element S(k) for a single s- 
wave entrance channel coupled to an isolated resonance 
scattering bound state, including a decay rate 7 to exit 
channels that lead to atom loss. The elastic scattering S- 
matrix element S(k) (equivalent to Eq. (3.13) of Ref. [371 ] J 
for an isolated decaying resonance is: 


S(k ) = 1 - 


iHT(k) 


E — E r , 


-i\h[^ + T{k)\ 


0 2iT) bg (k) 


• ( 21 ) 


where 


E res = h(co n + s n I — ui) = —HS (22) 

is the energy of the field-dressed molecular resonance 
level n. Its “bare” location at hui n = hoj a — E n is 
shifted by an intensity-dependent shift Hs n I , where E n is 
the binding energy with respect to the excited separated 
atom limit when 1 = 0. We define the molecular detun¬ 
ing S so it is negative for red detuning, in which case a 
resonance peak occurs when E = E res . 













We assume low power, in which case the shift varies 
linearly with intensity. The coefficient s n can be either 
positive or negative [371 IfiZl • where a positive value cor¬ 
responds to a shift of the resonance peak closer to the 
atomic line. The closed channel resonance bound state is 
coupled to the entrance channel by the stimulated decay 
rate, 

27T 

r(fc) = T MK pt |£}| 2 . (23) 


Here, | n) represents the excited bound state, which in 
general would be a mixture of the two |2) and |3) excited 
case (c) states in Table I. In practice |n) would be well- 
approximated by a single 0 U or 1 „ J = 1 vibrational 
state. The ground state scattering wavefunction \E) is 
assumed to be energy normalized [76], so that 




1/2 


sin(fci? + ?7bg) as R —> oo , (24) 


where the phase shift 7ft,g is related to the scattering 
length Ob g in the k —> 0 threshold limit as ? 7 bg = —ka^g. 

We emphasize that the form of the expression in 
Eq. (1211) is completely general for any isolated threshold 
resonance and applies equally well for MFRs and OFRs, 
if the terms are identified properly. The Fermi golden 
rule width T(fc) expresses the strength of the resonance 
pole term with a tunable denominator. When T(k) = 0, 
there is no Feshbach resonance, and we recover the stan¬ 
dard expression S(k) = e 2i7?bg for the uncoupled entrance 
channel. The expression in Eq. (© ensures that T(fc) 
follows the standard threshold law, and thus for an s- 
wave entrance channel, T(k) oc k. For non-decaying res¬ 
onances, 7 = 0 , and the imaginary term in the denomi¬ 
nator vanishes as k — > 0 . 


IV. COMPARISON BETWEEN OPTICAL 
FESHBACH RESONANCES AND MAGNETIC 
FESHBACH RESONANCES 

A. Isolated Resonance Theory 

The resonance length formalism summarized in Sec¬ 
tion II.A.3 of Ref. I] shows how to relate MFR and 
OFR resonance strengths and compare OFRs to MFRs 
in a unified approach. All we need to note is that in 
the case of an isolated MFR, the threshold width T(k) 
in Eq. m is given by an expression similar to Eq. (1231) , 
except that I4 P t needs to be replaced with an appropri¬ 
ate internal spin-dependent Hamiltonian mu. Further¬ 
more, E les would be replaced with a H-dependent tuning 
and shift mm, E res = - B c ) + E s hift, where 

B represents magnetic field, B c is the field where the 
bare resonance level crosses threshold, represents the 
difference between the sum of the magnetic moments of 
the two atoms and the magnetic moment of the bare res¬ 
onance level, and £’ s hift represents an energy-dependent 
shift term. 


The threshold law for s-wave collisions shows that as 
k —> 0 the quantity HT(k)/k (for either an OFR or the 
MFR analog) becomes a fc-independent constant having 
the units of length times energy. Thus, for either an MFR 
or an OFR, we can decompose T(k)/k into a product of 
a length factor L r and an energy E r , 


HT{k) 

2k 


L r E r , 


(25) 


Since only the L^E V product is significant, we are free to 
choose either the length L r or the energy E r factor to 
yield a convenient expression for the scattering length. 

In the case of non-decaying MFRs, it is conventional 
to choose L r = dbg- The E r factor is typically written 
as 8)jl A, thereby defining the magnetic “width” A of the 
MFR: 


= 0 be(< „a). (26) 

Here the subscript “MFR” indicates the type of reso¬ 
nance. When this form is substituted in Eq. ED, 7 is 
set equal to zero, and the k —> 0 limit is taken with the 
shift term in Refs. [H,[Z3], Eq. © reduces to the standard 
expression for an MFR, 


a (B) — abg Ubg g _ g , (27) 

where Bq = B c — £( s hift/<5/r is the pole position. In the 
case 7 7 ^ 0 , this procedure would give the complex scat¬ 
tering length for a decaying MFR )1,, [4l| . 

By analogy to MFRs, one can define an OFR resonance 
frequency width w by 

M- (28) 

Note that — ab g u; is positive definite. In the limit |<5| 

7 where we can ignore the decay of the resonance, the 
scattering length due to an OFR is 


w 

a(uj) = Ob g — Ubg 7 : pr ■ (29) 

W-K + s n I) 

The standard way to express the L r E r product in the 
case of an OFR is to define E r to be the known quantity 
tvy m and call the length parameter the “optical length” 

Zo P t [Hill, 


hr o F R(fc) _. , 

2jL — opt \k ,r )m) ■ (30) 

For large detunings \5\ 7 the scattering length is 


i(cj) — dbg + lopt~ 


r )m 

(oj n + S n I) 


(31) 


A similar resonance length parameter has been defined 
for a decaying MFR by Hutson by using the total 
decay width of the resonance for E r . 
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In the case of an OFR that decays only to the ground 
state, choosing E T = tv) m has the advantage of eliminat¬ 
ing the dipole strength from the expression for l op t- Using 
fry a = 327r 3 d 2 /3A 3 , the definition = 2q a , taking U opt 
in Eq. 0, and assuming the ^-independent molecular 
dipole moments of Eqs. m or m, we find 


, _ m 

° pt “ 2 k lm 


A l \{n\E)\ 2 
167rc k 


I fri 


(32) 


where / rot is equal to 1 for 0„ states and 2 for l u states 
(due to the different rotational Honl-London factors for 
parallel and perpendicular transitions). Consequently, 
L r E r is proportional to the product of a Franck-Condon 
factor and the square of the molecular electronic transi¬ 
tion dipole moment. 

Equation (13^1) shows that l opt varies linearly with PA 
laser power /. The only molecular physics parameter 
it depends on is the free-bound Franck-Condon factor 
|{n|E}| 2 , which varies linearly with k at small k. As an 
example, direct calculation of \(n\E)\ 2 /k for the J = 0 
n = -4 0„ level shows that this quantity decreases at 
a rate of 0.66% per /rK as E/ks ranges from 0 to 10 
/rK. Thus, l opt /I is only weakly dependent on collision 
energy in ultracold gases and may be approximated as a 
constant. Its weak variation with energy could be esti¬ 
mated from approximate theories based on the reflection 
approximation f37| or the stationary phase approxima¬ 
tion H. 

A useful way to compare the strengths of MFRs and 
OFRs is to use a dimensionless resonance “pole strength” 
parameter that applies to either case: s res = L v E T /aE, 
where a is the mean scattering length of the van der 
Waals potential [79| and E = h 2 /(2/j,a 2 ) is the corre¬ 
sponding energy. These are a = 71.8 a o and E/h = 
7.97 MHz for 88 Sr collisions. Chin et al. [lj used 
s res to characterize and classify MFRs according to 
whether s res > 1 (open channel dominated) or s res < 1 
(closed channel dominated), where the former tends to 
be “broad” and the latter “narrow.” Thus we have 


mfr _ Qbg ofr 

*res - ? *res 

CL jC/ 


Cpt kr/ rn 
a E 


(33) 


One obvious difference between MFRs and OFRs is that 
the strength of an OFR can be controlled by increasing 
the PA laser power to increase l opt , whereas the strength 
of a MFR is fixed. However, Z opt cannot be increased 
too much since the light scattering loss rate due to either 
atomic or molecular processes also increases with I 34). 

Experimentally useful MFRs tend to have a pole 
strength parameter between unity and 100 [lj. The width 
ratio frfm/E = 0.0019 is much less than unity for the 
narrow OFRs near the intercombination line of 88 Sr, so 
s®™ -C 1 unless it can be compensated by making Z opt /a 
very large compared to unity. Thus, 88 Sr OFRs tend 
to be weak, narrow, “closed channel dominated” reso¬ 
nances. An interesting comparison is with the experi¬ 
mentally useful broad open channel dominated MFR of 
two 85 Rb atoms at 155.2 G, for which s res = 28 [l| . This 


is a decaying MFR in an excited spin channel (38j, with 
a natural decay width of 7 /( 27 r) = 5.0 kHz due to spin 
relaxation of the “bare” closed channel state of the res¬ 
onance. The lifetime of 32 fis fic| of this spin channel 
is comparable to that of the excited Sr 3 Pi state. The 
major difference between the 85 Rb MFR and 88 Sr OFRs 
is the much smaller resonance strength s res of the lat¬ 
ter at intensities where the atomic light scattering is not 
harmful. 

In contrast to 88 Sr, OFRs for the species 172 Yb were 
found to have an Z Q pt on the order of 10 4 oo at I = 
1W/cm 2 for levels near 1 GHz atomic detuning [8fl| ■ This 
implies that broad open-channel-dominated OFRs with 
s res > 1 may be realizable with 172 Yb. It is not yet known 
whether OFRs might exist with s les > 1 for frequencies 
near the alkaline earth intercombination line in mixtures 
of alkaline earth species and alkali metal species. This is a 
subject for future experimental and theoretical research. 


B. Multiresonance Theory 

It is useful to compare isolated OFRs and MFRs since 
isolated resonance theory is widely utilized in both cases; 
however, both OFRs and MFRs exist as sets of reso¬ 
nances that interfere with one another, so it is instruc¬ 
tive to compare multiresonance treatments of the two 
effects. Although a thorough treatment of a multireso¬ 
nance OFR-MFR comparison could be the subject of an 
entire publication, in this section we provide an overview 
of such a comparison using results from multichannel 
quantum defect theory (MQDT). 

When resonance interference is considered, significant 
qualitative differences between OFRs and MFRs emerge. 
There are two sources of such differences. First, the 
molecular physics that determines the resonance strength 
is due to short range spin-dependent interactions for 
MFRs and long-range photoassociation for OFRs. Sec¬ 
ond, OFRs typically span many vibrational levels of the 
same closed channel molecular state whereas experimen¬ 
tally utilized MFRs are typically different spin compo¬ 
nents rather than a vibrational progression. 

Sets of overlapping MFRs are well-studied for differ¬ 
ent alkali metal species [U [78|, fsTI - lsil , and overlapping 
MFRs have recently been shown to be important for few- 
body physics [S5|. The interference of overlapping MFRs 
can be quantitatively explained by MQDT [78j, [8(1 [87], 
with which one can derive an S'-matrix that is a multires¬ 
onance generalization of the isolated resonance formula 
in Eq. (l2ll) . We introduce the MQDT theory here to 
highlight the differences between overlapping OFRs and 
MFRs. 

Considering one background channel and one closed 
channel, the background is characterized by the usual 
E-dependent phase shift ??b g (E) [78l I86l - l88l |. The closed 
channel c is characterized by a bound state phase func¬ 
tion v c (E — E c ). The energy E c of the separated atoms 
in the closed channel is modified with the “field tun- 
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ing,” which means varying the external magnetic field 
in the MFR case or the atomic detuning in the OFR 
case. Bound states of the closed channel exist where 
tan v c (E — E c ) = 0. Thus, the external field tuning 
moves the bound state spectrum relative to the back¬ 
ground channel E = 0 threshold, allowing bound states 
to be tuned across threshold. The coupling between the 
background and the closed channel is characterized by 
the dimensionless MQDT parameter s res , which may also 
depend on the external fields. 

If we follow Ref. (fU and express energies as e = E/E 
and k = fed, then the equivalent MQDT expression to 
Eq. (1211) for an MFR or an OFR can be written in a 
universal dimensionless form, 


Smqdt = <Sbg ( 1 — * 


2ks t . 


tan l/ c -f- C s hift T ^Cwidth 


where Sb g (e) = e 2l77bE ^-* is the background term, 


(34) 


e m (e) = 


dv c {e)/de 


(35) 


represents the mean spacing between different eigenener- 
gies, and e sh ift(e) and e w idth(e) = ^h'y / E + KS les represent 
the respective shift and decay parts of the complex energy 
of the interacting, decaying resonance. The shift term 
£shift(e) will scale as s res , and both of these quantities 
are only very weakly dependent on energy near thresh¬ 
old. Here we need to view the MQDT parameter s res as 
a continuous function of the external field that defines 
the Hamiltonian. 

Using Eq. (0, we can now describe some key differ¬ 
ences of OFRs and MFRs. These come from the variation 
with field strength of both the numerator and the denom¬ 
inator of Eq. (1341) . In the MFR case, it is an excellent 
approximation to take s res to be a constant, independent 
of B and E, since the interactions that determine this 
parameter are short range, where R -C a and the energy 
scale is large. Consequently, the matrix element that sets 
the magnitude of s res is independent of small field tun¬ 
ing. On the other hand for an OFR, s res scales linearly 
with laser power, and we must also think of s res (w) as 
being highly sensitive to field tuning, since the optical 
coupling is determined by the Condon point at very long 
range (on the order of a or larger). The Condon point 
varies rapidly with PA laser frequency, so the crossing 
structure of the field-dependent Hamiltonian varies with 
laser frequency in a major way, changing the response of 
the system to the optical field. Another way of think¬ 
ing about this variation is that for an isolated resonance 
s res is proportional to a Franck-Condon factor, which will 
vary rapidly from level to level in the closed channel, so 
that the general MQDT coupling parameter can not be 
taken as a field-tuning-independent parameter (89j . 

Secondly, note that the proper MQDT expression in 
the denominator of Eq. 0 that contains the effect of 
field tuning is the tan v c term, which vanishes at res¬ 
onance poles. To get the normal isolated resonance 


approximation near a tunable eigenenergy e res , as in 
Eq. (1211) , it is necessary to expand this function in a Tay¬ 
lor series as 

4- ( \ d " c ( e ^ c ) i / \ /oa\ 

tan u c (e - e c ) «--- £=£res (e - e res ) ■ (36) 

oe 

In the ultracold case, e tends to remain very close to 0 
but e res is varied by tuning the field, so the expansion 
in Eq. (1361) should be made near e = 0. This lineariz¬ 
ing approximation is normally quite good as long as the 
range of expansion e — e res remains small compared to 
the mean spacing e m to adjacent levels. This is normally 
the case for MFRs, where the widths of even broad reso¬ 
nances tend to be small compared to the distance to the 
next vibrational level in the same spin channel (90| . On 
the other hand, it is common to observe a series of vi¬ 
brational levels of the same electronic state in the OFR 
case. In this case, the tan v c must be left unexpanded 
and u c treated as a continuous function of field tuning if 
multiple resonances are present [11- 

Reference [78] showed how to extend the MQDT for¬ 
malism for Eq. (1341) to multiple spin channels with 
overlapping resonances. Generally, for the reasons dis¬ 
cussed above s res for each separate channel can be well- 
approximated as a E- and H-independent constant, and 
the linearizing approximation in Eq. (1361) is used for de¬ 
tunings spanning multiple spin channels. Then the gen¬ 
eralization of Eq. (1361) gives a sum of resonance terms 
similar to that in the pole term of Eq. (15S1) , where there is 
a global background scattering length ab g for all channels 
and the shift terms in each denominator depend on all the 
poles simultaneously. The formula can be transformed 
to a form where each narrow resonance can be viewed 
as an isolated resonance having a “local” (in field tun¬ 
ing) background modified from the global one by nearby 
interfering resonances. An extension to OFRs from dif¬ 
ferent electronic states may not be possible, because of 
the rapid variation of MQDT parameters with field tun¬ 
ing. Furthermore, the MFR theory should not be used for 
different members of the same vibrational series because 
of the inability to linearize the tan zz c function across two 
or more vibrational levels. It may be possible to de¬ 
velop some approximations appropriate to the OFR case, 
but meanwhile the isolated resonance approximation or 
CC calculations remain the best tools for understanding 
OFRs. 


V. ELASTIC AND INELASTIC COLLISIONS 
A. OFR isolated resonance formulas 

For detunings as large as hundreds of linewidths 
from the resonance center (but smaller than the separa¬ 
tion between resonances), the complex scattering length 
a(k,oj,I) and the elastic or inelastic collision rate co¬ 
efficients derived from the S'-matrix of Eq. m are in 
excellent agreement with the full CC calculations. In 
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this regime, Eq. m fully describes an isolated OFR as 
a function of energy, detuning, and intensity. 

The expression in Eq. ED gives the isolated resonance 
approximation to the entrance channel loss probability 
Pioss hr terms of only two parameters, Z 0 pt and the reso¬ 
nance position E les : 


Pioss = 1 — | S' (At) | 2 


where 


D 2 + \(l + 2kl eS ) 2 ’ 


^eff — ^opt(bm/b)? 

D = (E — E res )/h'y. 


(37) 


(38) 

(39) 


Here the effective optical length l e g determines the 
resonance coupling strength for the general case when 
7 > 7m- -Pioss determines the inelastic cross section in 
Eq. © and thus the loss rate coefficient K in in Eq. @. 
Note that for 2 kl e g 1, Eq. [37] describes power broad¬ 
ening. 

Figures ED a) and 0b) show the behavior of Pi OS s as a 
function of detuning for different intensities and collision 
energies. In Fig. ©a), P\ oss is not saturated at low colli¬ 
sion energy. As collision energy increases, Poss broadens, 
and its peak value (when E = E les ) approaches its upper 
bound of unity. This only occurs for red molecular detun¬ 
ings. Fig. 0b) shows similar broadening in the inelastic 
rate coefficient K in oc Pi oss as intensity is increased and 
collision energy is kept low. For large intensities, Ki n 
saturates at its value given by the unitarity limit. Note 
that according to Eq. [22l the intensity-dependent fre¬ 
quency shift of the resonance is accounted for in 6. For 
temperatures in the nK range, thermal averaging of the 
line shape is essential to compare with experiment. As 
is well-known [l9i | . PA lines have a pronounced red-blue 
asymmetry when ksT is larger than the natural width 
of the PA line. 

The isolated resonance S matrix can be used to derive 
the complex ^-dependent scattering length a(k). Com¬ 
bining Eqs. ET1) and JTJ, 


a(k) 


Qhg(k') T 


□gFl ( X + k 2 a he {k) 2 ) 

E - E res - ka hg (k)^ + ’ 


(40) 


where «bg (k) is found by using the background Sb g (k) = 
g 2 *r?bg j n gq_ 0 - Notice that this expression does 
not contain power broadening, which enters the elastic 
and inelastic cross sections through the f(k) factor in 
Eqs. 0-@. 

In the limit that fc|ab g | ■C 1 (valid for a Bose-Einstein 
condensate), we obtain 


a a hg +l e s 62 + ^ 2 j 4 2^ eS S 2 + y 2 /4 ’ ^ 

where we have taken a^ g = Ob g to be real. Fig. [5] illus¬ 
trates the real and imaginary parts of a = a — ib calcu¬ 
lated at E/ks = 1 nK. The results from Eq. (Ell) are in 




FIG. 4: (Color online) (a) Coupled-channel calculated in¬ 
elastic loss probability 1 — |5(fc )| 2 versus atomic detuning 
v — vq for different values of collision energy E/ks- We have 
used 1 = 5 W/cm 2 and the 0 U J = 1 n = —4 feature, 
where E n /h = —1084.0763 MHz. The black dots indicate 
the peak values calculated with l op t =161.5 ao, 7 = 7 m , and 
v — i/o = ( E n — E)/h (where D vanishes in Eq. (1371) '). Since 
1 — |S(fc )| 2 calculated from the isolated resonance formula in 
Eq. (1371) is indistinguishable from the CC calculation on this 
graph, only the peak comparisons are shown by the black dots 
on the figure. Panel (b) Coupled channels calculated inelastic 
loss rate coefficient Ki n versus molecular detuning. Here we 
plot different values of PA laser intensity I. We use a collision 
energy E/kB = 1 nK and the 0 U J = 1 n = —4 feature near 
v — vq = —1084 MHz. The black dots show the predictions 
of the analytic formula in Eq. (ED- The line wings beyond 
around |<5/7m| = 6 scale linearly with power. The black dot¬ 
ted line represents the unitarity limit where Kin saturates. 


excellent agreement with the CC calculations. The peak 
values of b = 2 Z opt at S = 0 and of a = ±/ op t at S = ± 7/2 
are also plotted in Fig. [5] 
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FIG. 5: (Color online) Coupled channel calculation of a and 
b versus molecular detuning in line width units 5/7 m, with 
7 = 7 m . We use the 0„ J = 1 n = —4 feature and 10W/cm 2 
PA laser intensity. The dots show the analytic predictions at 
the peak extrema using Eq. (S3 with lopt — 323.2 no- 



FIG. 6 : (Color online) Trap loss, given by Eq. 1451 as a func¬ 
tion of detuning for various cloud temperatures. Here 7 = 7 m , 
Z eff = lOOao, tpA = 200ms, and we take the limit of r > tpA- 
We also assume a spherical trapped cloud with No = 6 x 10 4 
atoms and a 20 /tm r.m.s. radius. Note the broadening to¬ 
ward red detunings, which is a result of momentum averaging 


B. Photoassociation 

The isolated resonance approximation also describes 
photoassociation [92j , the process by which two ground 
state atoms and a photon combine to form an electron¬ 
ically excited molecule [lt|. This phenomenon can be 
used to measure l e g, which characterizes the strength 
of an OFR. This strength can be inferred from mea¬ 
surements of the trapped atom loss (into untrappable 
molecules) that results from driving a photoassociation 
resonance. Using Eqs. m and ©, the inelastic rate 
coefficient, which describes molecule formation, is 


Kj/n (4, kJ) 


471 -h _ 7 2 leS _ 

m {5 + E/hy + ^(i + 2ki eS y 


(42) 


For a trapped ultracold thermal gas, one must intro¬ 
duce Boltzmann averages into the theory. To this end, 
we approximate that the PA laser interacts with an entire 
velocity class at each point in space within the trap. This 
approximation, which is good for large densities, means 
that one must momentum average K-i n . Furthermore, we 
use the fact that photoassociation is a short-range phe¬ 
nomenon in the isolated resonance approximation j37| ; 
therefore, we do not perform a Boltzmann spatial aver¬ 
age in this treatment. 

The quantity Ki ni which is the momentum-averaged 
K in , is given by 


K — 

1 m — 


L m q a 

7r d P5 


V^ k ti 


"d 3 p 1 Jd 3 p 2 e -^+^)/Po K in (6, l e s, k) 

/»00 

/ dkk 2 e- k2 / k * h K m (6,l e s ,k), (43) 
Jo 


where po = \j2mkBT , k t n = yJIpksT / H, p\ and p 2 are 
the momenta of the two collision partners, and \p\— p 2 \ = 
Hk. 

In order to use K in to describe trap loss due to 
photoassociation, we must understand what fraction of 
molecules is ejected from the trap. As mentioned in 
Section III B1 we approximate that every photoassociated 
molecule is lost to the trap. We have numerically studied 
this approximation, finding that it is good to 1% for all 
molecular states except for the least-bound 0„ level. In 
this case, the evolution of the in-trap atomic density is 


p=-2K in ^-- P -, (44) 

Z T 

where p is the atomic density and r is the one-body life¬ 
time (due to loss mechanisms such as background gas 
molecules and atomic light scatter). Here the p 2 term 
arises from the number of pairs in an IV-particlc sam¬ 
ple, N(N — l)/2 ~ N 2 / 2. The signal in photoassociation 
experiments is typically the atom number N after the 
application of the PA laser (H, [63j, given by integrating 
the solution to Eq. (l44l) over space, 


N(6, l eS ) 



Po(r) e tpA / T 

1 T A i n ((5, Z e ff) po(t') 7"(1 




(45) 

Here tpA is the pulse duration of the PA laser and po{r) 
is the in-trap density before the PA laser is applied. Fig. 
© depicts N(5,l e ff) for different temperatures and de¬ 
tunings. The density po(r) can be determined by fitting 
experimental in-trap absorption images to a 3D Gaussian 
distribution (Ref. |34(, supplementary online material). 

Unless a magic wavelength trap is employed, the ac 
Stark shift from optical traps causes a position-dependent 
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laser cooled 88 Sr) is greater than the PA laser photon re¬ 
coil temperature, Doppler broadening is negligible com¬ 
pared to the momentum broadening shown in Fig. [ 6 ] We 
numerically checked whether Doppler effects are signifi¬ 
cant for our analysis, and we find that Doppler broad¬ 
ening can only be neglected in the fits of Fig. [7] due 
to the presence of both a large momentum broadening 
and an appreciable Stark shift broadening toward blue 
detunings. 

Treating collisions in this manner breaks down when 
elastic processes become important. The elastic-to- 
inelastic collision ratio is approximately 

K el (k)/K m (k) ~ 2 kl eS = 2 k th l eS , (47) 


v - v 0 (MHz) 


FIG. 7: (Color online) Atom number data from Ref. [34] as a 
function of the PA laser detuning from the atomic resonance. 
The false color pictures above the plot are the measured opti¬ 
cal depths corresponding to the data points directly below the 
centers of the pictures. These measurements were performed 
in a non-magic wavelength trap, which caused broadening to¬ 
ward blue detunings since LOstark{0) = 2-zr x 327 kHz. The 
solid line is a fit using Eq. (14511 with LUstarkir) included in 
Ki n . The quantities T and l e fj were free parameters in the 
fit. 


frequency shift of the atomic resonance [HI, [63j. For pho¬ 
toassociation experiments, this effect results in a broad¬ 
ening of the lineshape feature in the signal N. To model 
this broadening, one must understand both the atomic 
response to the optical trap and the intensity profile 
I trap (p) of the trap laser (Ref. j34[, supplementary online 
material). In the case of 88 Sr, the 1 Sq and 3 Pi polariz¬ 
abilities are known well enough to calculate the differ¬ 
ential polarizability, 013p 1 (u/ t rap) ~ ^Soi^trap) > to better 
than 10% for typical trap laser wavelengths (such as 1064 
nm). 

One can model Itrapif) with the Gaussian beam equa¬ 
tion, using parametric heating measurements to obtain 
the beam waists {93|. This broadening can be included 
in the photoassociation signal by adding a Stark shift, 


^ Stark 


(r) = - 


2 HenC 


[o'3p 1 (^UJtrap) ^So (.^trap)\ ktrapir) > 


(46) 

to <5 in Eq. [42] and then carrying this term through to 
Eq. (l45j) . With the trap ac Stark shift accounted for, we 
can fit experimental photoassociation data (Fig. Q. For 
these fits it is necessary to approximate the integrals in 
Eqs. m and (H51) as sums (Ref. 0, supplementary 
online material). 

A full treatment of the photoassociation lineshape 
would include Doppler broadening. However, accord¬ 
ing to theoretical studies of narrow-line photoassocia¬ 
tion [92l |. since a T of a few fj,K (typical for narrow-line 


where the overline denotes thermal averaging. In deriv¬ 
ing this formula, we made the approximation e 2lVb3 ~ 1 , 
which is acceptable for the above estimate since ctbg is 
only -1.4 ao- Therefore, when l e g ~ l/2kt.h, elastic colli¬ 
sions must be treated. 


C. Elastic Collisions 

If photoassociation can be minimized and elastic col¬ 
lision rate modifications can be made large, the OFR 
effect could be used to manipulate atomic interactions 
with relatively little particle loss. The usefulness of such 
manipulations is evident from experiments based on the 
MFR effect, which has proved to be a very fruitful tech¬ 
nique that is central to numerous experiments [H. To 
access this regime in a quantum degenerate gas (for 
which k —> 0 ), the ratio of the optically modified elas¬ 
tic scattering length to the inelastic scattering length, 
[a( 0 ) — abg\ / 6 ( 0 ) = 26 / 7 , must be much greater than 
unity. For a thermal gas, the rate coefficients determine 
the relevant limit, (I\ e i(k)/I\i n (k)) 1 , which implies 
via Eqn[47|that l e s 5S> 1/2 k t h- 

We estimate from the latter condition that a thermal 
88 Sr gas at T = 3 fj,K (typical of narrow-line laser cooling) 
will require l e s to be much greater than 400 ao- Using 
the isolated resonance approximation, large changes in 
elastic scattering were predicted to arise from an OFR 
based on the n = — 1 vibrational state [30]. This pre¬ 
diction required high PA laser intensity and very large 
molecular detunings from the n = — 1 state. However, as 
our CC theory has shown, these conditions will not yield 
an effect comparable to a lossless MFR because the req¬ 
uisite detunings are larger than the separation between 
resonances. Instead, the OFR physics is determined by 
the nearest resonance to a given detuning ( Section |II Dl) . 

Experimentally, elastic collisions can be studied using 
cross-dimensional thermalization. For instance, in Ref. 
f94j, a trapped atomic gas was prepared in a nonequilib¬ 
rium state using parametric heating in ID of the trap. 
Due to elastic collisions, the authors observed thermal¬ 
ization of the non-equilibrium gas. The thermalization 
time of the gas can be calculated from a simple treatment 
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v-vo (MHz) 

FIG. 8: (Color online) Temperature data from Ref. [m| for 
horizontal (H) and vertical (V) trap eigenaxes. The solid lines 
are the results of a Monte Carlo simulation including elastic 
and inelastic collisions. 

based on Enskog’s equation or a full molecular dynamics 
simulation [95[. These treatments show that, on average, 
each particle participates in about three elastic collisions 
events during the 1 /e thermalization time. 

Elastic collisions induced by OFRs were experimen¬ 
tally studied in Ref. [34| using cross-dimensional thermal¬ 
ization. In this work, an OFR was accessed in a trapped 
88 Sr gas prepared in a non-tlrermal state. In the absence 
of OFR-induced collisions, this gas would not thermal- 
ize over experimental time scales because of the negli¬ 
gible dbg = —1.4ao in 88 Sr. With a PA laser applied, 
clear temperature changes were observed as a function of 
v — v q. To understand whether these observations were 
caused by elastic collisions, we apply our theory to the 
experimental data. 

Since the atomic samples in this measurement never 
reached thermal equilibrium over the time scale of the 
experiment, it is not possible to study the resulting data 
analytically. Instead, a time-dependent simulation of the 
phase-space density is necessary to understand the data 
quantitatively. Our analysis uses a numerically efficient 
method due to Bird [96] for simulating collisions between 
thousands of particles. The method discretizes the trap 
volume into small “collision volumes” containing much 
less than one particle on average. If there is more than 
one particle in a collision volume, the probabilities P e i 
and Pi n of elastic and inelastic collisions between these 
atoms is calculated as P e i/m = \vi — ubl <J e i/intstep/V■ 
Here V\ and V 2 are the velocities of the two atoms, t step is 
the time step of the simulation, V is the collision volume, 
and the cross sections are given by 

8tt [ksl + a bg (S + §)] + a b g 7 2 /4 //|cA 

1 + (S + E/h) 2 + \ (1 + 2 kleff) 2 

_ 47T leg _7^_ 

k {5 + E/h) 2 + ^(l+2kl eS ) 2 ' 


If an inelastic collision occurs, both particles are re¬ 
moved from the simulation. For an elastic collision event, 
the particles’ velocity vectors are rotated using a ran¬ 
dom rotation matrix. Between each of these collision 
steps, the particles evolve in the trap potential using an 
embedded Runge-Kutta method. We have checked this 
simulation against known results for cross-dimensional 
thermalization in harmonic traps and have also confirmed 
that, with elastic collisions removed, the simulation re¬ 
produces the results of photoassociation theory of the 
previous section [97]. 

Fig. [ 8 ] de picts temperature data from the experiments 
of Ref. [3J] as well as our simulation results. Our sim¬ 
ulation tells us that the temperature peaks for certain 
detunings because the PA laser is driving photoassocia- 
tive loss of the coldest atoms, resulting in antievaporative 
heating. We also find that without including elastic col¬ 
lisions in our simulation, the simulation does not predict 
the dip in the horizontal temperature apparent in the 
data. Therefore, we conclude that this temperature dip 
indicates partial thermalization of the gas. The fact that 
antievaporation and thermalization have different detun¬ 
ing dependence arises because the elastic and inelastic 
collision rates average differently over the collision mo¬ 
mentum k. 

The interplay between elastic and inelastic processes is 
sensitive to the value of 7 used in the simulation. The 
simulation only agrees with the experimental data if we 
set 7 = 27t x 40 kHz = 2.7^ m . This leads us to conclude 
that the OFR effect in 88 Sr is broadened beyond the nat¬ 
ural spontaneous decay of the 88 Sr 2 molecules. Extra 
broadening has also been seen in other 88 Sr mil and 
Rb [3l| OFR experiments. 


VI. SUMMARY AND CONCLUSION 


We have developed CC and isolated resonance theo¬ 
ries of OFRs. The CC theory predicts resonance inter¬ 
ference for detunings between OFRs, causing the OFR 
effect to vanish between resonances. We conclude that 
OFR experiments have a “nearest resonance” constraint, 
meaning that the OFR effect is dictated by the nearest 
resonance to a given detuning. The isolated resonance 
theory agrees with the more complete CC theory for de¬ 
tunings near a molecular resonance. In this regime, it 
is possible to use the simpler isolated resonance theory 
to model photoassociation and OFR measurements and 
fit the data from these experiments. Such models have 
shown a broadening beyond the expected linewidth of the 
molecular state accessed by an OFR. 


( 49 ) 
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